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ABSTRACT 

We present a new algorithm designed to improve the signal to noise ratio (SNR) 
of point and extended source detections in direct imaging data. The novel part of 
our method is that it finds the linear combination of the science images that best 
match counterpart images with signal removed from suspected source regions. The 
algorithm, based on the Locally Optimized Combination of Images (LOCI) method, is 
called Matched LOCI or MLOCI. We show using data obtained with the Gemini Planet 
Imager (GPI) and Near-Infrared Coronagraphic Imager (NICI) that the new algorithm 
can improve the SNR of point source detections by 30-400% over past methods. We also 
find no increase in false detections rates. No prior knowledge of candidate companion 
locations is required to use MLOCI. While non-blind applications may yield linear 
combinations of science images which seem to increase the SNR of true sources by 
a factor >2, they can also yield false detections at high rates. This is a potential 
pitfall when trying to confirm marginal detections or to re-detect point sources found 
in previous epochs. Our findings are relevant to any method where the coefficients of 
the linear combination are considered tunable, e.g. LOCI and Principal Component 
Analysis (PCA). Thus we recommend that false detection rates be analyzed when using 
these techniques. 
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INTRODUCTION 


While more than 1700 planets have been detected by transit and radial velocity techniques, 
only about a dozen planetary-mass objects (< 13 Mj up ) have been directly im aged ar ound stars to 
dat^} However, this number is expected to increase significantly (e.g. [McBride et al. 2011) since 


first light has just been obtained with the new extreme adaptive optics instruments, Gemini Planet 


Imager(GPI; Macintosh et al. 2008) and SPHERE (Beuzit et al. 2010). A severe challenge in AO 
imaging arises from quasi-static speckles, which can emulate astronomical point sources for hour 
long timescales (e.g. Marois et al. 2005). Besides high-order AO corrections and coronagraphy, the 
most powerful advances to overcome this sensitivity barrier are (1) angular differential imaging 
(ADI; Liu 2004| |Marois et al. 2005) which decouples the sky rotation of the planet from the 
speckles, and (2) spectral difference imaging (SDI; |Racine et al.|[T999| |Biller et al.| |200T|) and 
spectral deconvolution (Sparks &; Ford 2002| |Thatte et al. 2007) which decouple the star and 
planet spectrally, taking advantage of the fact that speckles move out radially with wavelength 
and that planet and stellar spectra differ. To obtain the maximum SNR for true companions, the 
speckle structure, which is correlated over time and wavelength, must be subtracted. Thus, some 
science images must take negative weight in the linear combination to create the final reduced image. 
Tremendous effort has been dedicated to finding algorithms that optimize detections by varying the 


weights of science images. The most prominent of these have been the LOCI algorithm (Lafreniere 
et al.||2007 ), and the PCA based algorithms, KLIP ( Sonmmer et al.||2012 ) and PynPoint ( Amara fe 
Quanz||2012 ), which have focussed on minimizing speckle residuals. Subsequent improvements have 
concentrated on optimizing SNR of point source recovery, e.g. damped LOCI, ALOCI and TLOCI 
( Pneyo et al.||2012 ; Currie et al.||2013 ; Marois et al.||2014 ). 


Our own experiments with data from the Gemini NICI Planet-Finding Campaign (Liu et al. 


2010) based on blind recovery of simulated companions showed no measurable improvement in 
achieved contrast using LOCI compared to conventional ADI (Wahhaj et al.| 2013a). We found that 
the main problem with LOCI was that reductions in noise around detections were also accompanied 
by reductions in the signal. Here, we present a LOCI variant that significantly improves SNR. We 
demonstrate the effectiveness of the algorithm using NICI Campaign datasets and the early science 
GPI dataset on HD 135344 B, a transition disk with strong evidence for ongoing planet formation. 
The science results from the GPI observations will be published in Menard et al. (2015). 


2. MATCHED LOCI 

The ADI reduction method is based on the principle that when the image rotator is turned 
off for an altitude-azimuth telescope, at the Cassegrain focus, the sky rotates with respect to 
the detector, while the quasi-static speckle pattern caused by the telescope and instrument optics 
do not. At the Nasmyth focus, one achieves the same effect in an observing mode called pupil¬ 
tracking. Astrophysical sources are thus separated from the speckle pattern, which is isolated by 
median combining a series of speckle-aligned images to create a reference PSF. The reference PSF 


1 http://exoplanet.eu/catalog/ 
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is then subtracted from each science image. The difference images are then de-rotated to align the 
sky and stacked again to create the final reduced image. The details of our own implementation 
of regular ADI can be found in Wahhaj et al. (2013a). The LOCI algorithm differs in the PSF 
subtraction step where it subtracts the best linear combination of speckle-aligned images such that 
the RMS in some region of the Field of View (FoV) is minimized. 

Below we list the steps of our new algorithm which we call MLOCI or Matched LOCI: 


1. A basic ADI reduction is initially applied to identify all candidate point sources with SNR>2. 

2. For each candidate, we define an annulus with inner and outer radii bounding the candidate. 
In this annulus, a sector subtending 40 degrees is defined with the candidate azimuthally 
centered and radially one PSF FWHM farther out than the inner edge. The rest of the 
annulus, subtending 320 degrees, is defined as the reference sector. 

3. For each science image, we create a counterpart where a point source is subtracted from two 
(or three) sky locations in the reference sector (see Figure [I]). The point sources which we 
subtract are simulated to undergo ADI self-subtraction, but for 30% less rotation as would 
occur in the ADI data set. However, the sources are not allowed to shrink azimuthally to 
widths <0.7xFWHM. The radial width of the sectors, the number of point sources subtracted 
and their intensity, we have found by experiments which we describe in §4. 

4. Next, LOCI is used to linearly combine the original science images to find the best match to 
each counterpart image to make counterpart PSFs. The science images should be speckle- 
aligned at this stage. When there are several spectral channels, the images need to be spec¬ 
trally deconvolved as described in §3. The matching is done only over the reference sector, 
not the target sector which includes the candidate location. 

5. The science images and their counterpart PSFs are then derotated to align North Up. 

6. An image with two (or three) positive sources at the sky locations of the previous negative 
sources (see Step 3) is made. The background of this image is zero i.e. it is noiseless (Figure [2j 
top row, middle) . 

7. Then LOCI is used to match the positive source image by linearly combining the derotated 
science images and counterpart PSFs (see Figure [2]). When there are hundreds of images 
(too time-comsuming for LOCI), we process subsets of the images and median combine the 
results. 

8. The sky locations of the injected sources are changed to positions not over-lapping with the 
previous ones. The steps above are repeated to obtain 5-10 different reductions. Since the 
injected sources are only needed to compute the coefficients of the linear combination, for the 
final combinations we use the original images without the source injections. Thus, these 5-10 
reductions only contain real astronomical sources and not the injected ones. 

9. The reductions above are then median combined to create the final reduced image for a 
particular annular region. 
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10. For X candidates, X annular reductions are made, the widths of which are at least as large as 
those of the matching annuli. The outer radius of these annuli are extended up to the matching 
annulus (inner radius) of the candidate next in radial separation. When the annular reductions 
overlap due to densely packed candidates, the overlapping pixels are median combined. 

It is important to remember that MLOCI only uses injected point sources to solve for the linear 
coefficients for image combination. Point sources are not injected at the locations of candidate 
companions but in the reference sectors which exclude the candidates. 

We also present here an implementation of MLOCI for extended sources, which one could 
apply if there is any evidence for extended emission in a basic ADI reduction. However, in this 
implementation the occurence of false positives becomes a concern. Therefore, an MLOCI detection 
of a data set devoid of signal should be performed to estimate the likelihood of false positives. 
Moreover, as a test, one can attempt to recover injected extended sources in empty data sets using 
MLOCI. This method is most promising when a library of PSFs from a different star are available 
as references, thus allowing one to avoid self-subtraction of the extended source. However, such 
detailed analyses are defered to a later publication. 

Below we list the steps of MLOCI for extended source detection: 

1. A preliminary reduction (e.g. basic ADI) is done first to identify all regions likely containing 
signal (e.g. SNR > 2). These regions are defined as the target, while the regions which are 
likely devoid of signal are defined as background (e.g. SNR < 1). These regions together 
make the reference region and will be matched by MLOCI. Other non-reference regions in the 
image will not be matched. Although, this step depends on prior estimates of the emmisive 
and non-emmisive regions around the target, we show that the final reduction is not very 
likely to yield false disk detections (see §4.3). 

2. For each science image, we create a counterpart where the local RMS is subtracted from the 
target sky regions. 

3. Then LOCI is used to linearly combine the original science images to find the best match 
to each counterpart image to make counterpart PSFs. As in the point-source reduction, the 
science images need to be speckle-aligned at this stage. The matching is done only over the 
reference region. 

4. The science images and their counterpart PSFs are then derotated to align North Up. 

5. The preliminary reduction is then set to zero at the background sky regions to make the final 
counterpart. 

6. Then LOCI is used to match the final counterpart by linearly combining the derotated science 
images and counterpart PSFs. When there hundred of images (too many for LOCI), we 
process subsets of the images and median combine the results. 
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3. OBSERVATIONS AND DATA REDUCTION 


We observed HD 135344 B in the J-band on April 21st, 2014 UT as part of GPI Early Science 
(Program ID = GS-2014A-SV-402). HD 135344 B is a prime example of a transition disk suspected 
to harbor a substellar companion (a brown dwarf or a recently formed giant planet). Submillime¬ 
ter images show a dust-depleted cavity with a radius of ^40 AU in a nearly face-on disk, while 
polarimetric differential imaging in the near-IR reveals two remarkable spiral arms extending from 
^28 AU (inside the submillimeter cavity) to ^130 AU (Garufi et al. 2013). 


We obtained images from the integral field spectrograph (IFS) in 37 spectral channels from 
A=l.l-1.35 /im using GPFs unprecedented adaptive optics capabilities. The star was placed behind 
a coronagraphic mask of radius 92 mas and thus a properly detected point source has projected 
separation, p>112 mas (since J-band resolution is ^40 mas ^|. The IFS has a plate scale of 
14.3 mas/pixel and a 2.8 // x2.8 // field of view (FoV). The telescope rotator was turned off allowing 
the sky to rotate on the detector through 35°, so that the speckle pattern produced by the telescope 
and instrument optics would be decoupled from real astronomical sources. This ensures a total 
azimuthal sky motion of 1.7xFWHM at 112 mas (and 2xFWHM at 130 mas). We obtained 39 
sixty-second exposures. A total of 50 minutes of telescope time excluding acquisition was used. 
A speckle at 112 mas should move outward by only (^p — l)xll2 mas=25.5 mas (<lxFHWM) 
across images obtained over A =1.1-1.35 /im and thus we do not expect much contrast gain from 
spectral deconvolution at this separation. 


We used the GPI pipeline (Maire et al. 2010) to produce wavelength-calibrated spectral cubes 
with bad-pixel removal, de-striping, flat-fielding, non-linearity, and persistence corrections. The 
centroid of each of the 37x39=1443 images are found by averaging the coordinates of the peaks of 
the four GPI satellites spots [^J The ADI datasets from each spectral channel are reduced separately 
using the pipeline described in Wahhaj et al. (2013a). However, we do not spatially filter the science 
images. Also, the reference PSFs are not translated to optimize the fit to the science images before 
subtraction, but only scaled in intensity to minimize the RMS between p=70-210 mas. The 37 
ADI reductions from each of the spectral channels are then median combined. 


For the MLOCI reduction, we create counterpart PSFs for each of the 1443 science images. 
Before combining science images to create the counterpart PSFs, they have to be spectrally de¬ 
convolved. That is, to match an image in the ith channel, the science images are magnified by a 
factor A^/Aj, where A i and A j are the central wavelengths of the zth and j th channels, respectively. 
The science images and their PSF counterparts images are then derotated to align North Up. 
We then use LOCI again to match the adjusted preliminary reduction, by linearly combining the 
science images and the counterpart PSFs. We have two counterpart PSFs for each science image: 
one made from images in the same spectral channel, and the other made from images with the 
same orientation on detector. In the last step, the 1443x3 science and counterpart PSF images 
can be combined using LOCI to match the final counterpart (See §2). However, this is often too 
many images for LOCI to combine. So instead, we sequentially combine all images (science and 


2 http: / / www.gemini. edu / sciops / instruments / gpi / instrument-performance? q=node /12172 

3 http: / / www.gemini.edu/sciops / instruments / gpi / public-dat a/pub lic-data-readme 
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two counterpart PSFs) in the j th channel and jth orientation with (37+39) x3=228 images in each 
combination. See Figure [3] for the reduction of the GPI data. 


We also demonstrate the effectiveness of the MLOCI algorithm using 10 methane band 


(CH 4 S, A =1.578 /im) datasets from the Gemini NICI Planet-Finding Campaign (2008-2012; Liu 


et al.pOlO ; Wahhaj et al.||20l3b ; Nielsen et al.||2013 ; Biller et ah1|2013 ). The filter was designed to 


help detect methane which is usually found in the atmospheres of cool substellar companions with 
surface temperatures <1400 K. In this paper, the simultaneously exposed images in filter CH 4 L 
are not used for spectral differencing as we want to keep the analysis simple. However, since we use 
a single narrow-band (4% width) filter, we reach shallower contrasts than in the Campaign. The 
Campaign targets used were HD 107146, HD 25457, HD 53143 , GJ 388, HD 31295, HIP 25486, 
HD 92945, HD 21997, HD 110058 and UY Pic, all observed between January and March of 2009. 
NICI was installed at the 8.1 m Gemini South Telescope and produced AO corrected images over a 
18.4" xl8.4" FoV detector. The NICI plate scale was 17.96 mas/pixel. The Campaign observations 
used a translucent coronagraphic mask (central attenuation factor=358; Wahhaj et al. 2011) with a 
half-power radius of 0.32". Median contrast achieved at p— 0".5 was 12.6 mag ( Wahhaj et al.|2013a ). 
We already demonstrated in Wahhaj et al. (2013a) that the LOCI pipeline does not perform better 
than regular ADI for NICI data. 


We compare MLOCI reductions of the NICI datasets to both PCA reductions using the KLIP 
algorithm (see Soummer et al. 2012) and regular ADI reductions using the implementation described 
above. The PCA algorithm creates an orthogonal basis of images which are then linearly combined 
to emulate the supplied PSF images. Moreover, the basis are ordered such that the first basis 
or mode is the best general match to the majority of the images, while each subsequent mode 
represents perturbations of decreasing importance. Thus fewer modes can be used in the hopes of 
reducing noise without removing too much signal, and improve over LOCI. The PSF subtraction 
can also be optimized for sub-regions of a science image. For our PCA reductions, we obtained the 
best results (best recovery SNR for injected companions) for number of KLIP modes = 8, and PSF 
subtraction optimized for annuli of width 20 pixels. The inner radius of the innermost annulus was 
set to 25 pixels. 


4. RESULTS 

4.1. Recovery of Simulated Companions 

To tune MLOCI and measure its optimum performance, we insert simulated companions into 
our 10 NICI datasets and recover them using MLOCI. We place the companions at 0".5-l".0 
separations from the primary (where sensitivity is speckle-limited) with contrasts just below the 50- 
detection limit for regular ADI reduction. For each data set between 0".5-l".0, we have roughly 800 
independent sky locations at which to do this test, given that the FWHM of the NICI iL-band PSF 
is ^55 mas. As described in §2, injected test companions are bounded by target sectors, while two 
or three point sources are inserted into the corresponding reference sectors (annlus — target sector) 
to aid MLOCI in preserving point source signal. Our aim is to compare the SNR of simulated 
companions as recovered by MLOCI, to the SNR as in recoveries by ADI and PCA reductions. We 
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measure MLOCI’s false positive rate by repeating the reductions without simulated companions at 
the test locations. 


First, we find the best reduction parameters for MLOCI given different amount of sky rotations 
in an ADI dataset. We multiply the true position angles (PA) of the FoV by a negative number 
to misalign the images as done in Wahhaj et al. (2013a) so that any real astronomical sources will 
be mostly removed in a median combination of the images. To simulate different amounts of sky 
rotation, we choose 5 values of this negative number so that the total sky motions at the separations 
of interest, which we call cp are 0.35, 1.25, 2.15, 3.1 or 4.0 xFWHM. We had to use companions 2 
and 0.5 mag brighter for a — 0.35 and 1.25 xFWHM respectively, because of reduced sensitivity 
in these cases due to large self-subtraction factors. We also vary W, the width of the reference 
sectors, over 4, 8, 16, 32 and 64 pixels and P the pull over 1,2,3,4 and 5. Here, P is the brightness 
relative to the local RMS, of the point sources that are subtracted from the reference sectors (see 
§ 2). Reductions of the UY Pic dataset by varying the parameters over this grid of values show 
that optimum settings are near W =16, P =3, while two simulated companions in the reference 
sector to assist MLOCI are adequate. 


Using the above parameters, we show two MLOCI reductions (for a =0.35 and 1.25) of the 
UY Pic (K0V, H— 5.9 mag, age ~ 70 Myr; Lopez-Santiago et al. 2006) data set along with ADI 
and PCA reductions (see Figure [5]). In these reductions, simulated companions were placed at 
separations 0''.5, 0A6, 0''.7 and 0''.8 and at PAs 360, 270, 180 and 90° just below the detection 
limits of an initial ADI reduction. The SNR improvement factors achieved for MLOCI over ADI 
and PCA range from 1.3 to 2.3. 


In the SNR calculations, for the noise we use the RMS in an annulus of width 2xFWHM 
centered on the test locations but exclude circular regions of diameter 1 xFWHM around them. 
The signal is just the intensity at the test locations minus the median intensity in the noise region. 
To decide whether a source in the reduced images is a detection based on shape, we calculate the 
fractional reduction in the RMS in a box of size 2xFWHM+l pixels when an optimally scaled 
Gaussian PSF is subtracted from the source. Here FWHM is estimated for the source, first. When 
the fractional reduction is >0.3 and the SNR is 5, we consider the source a detection, a criteria 
which recovers a 5 a signal >99% of the time with a negligible false detection rate (FDR; sometimes 
termed false discovery rate) in our tests. We use this criteria over that in Wahhaj et al. (2013a) 
because it has a clearer interpretation. 

In the example MLOCI reductions of Figure [5j no other sources besides the ones which were 
injected into the data were recovered according to the detections criteria specified above. Although 
MLOCI can improve SNR by a factor >2, we need to check that the FDR is not elevated compared 
to ADI. The FDR is the fraction of detections that are incorrect, i.e., they do not result from real 
or injected signal. It is important that this number is low, so that negligible time is wasted in 
follow-up observations trying to confirm false detections. The FDR is also different from FPR, the 
false positive fraction (see Figure [4]). Next, we discuss FDR for MLOCI reductions. 


Wahhaj et al. (2013a 








4.2. Contrast, SNR, Completeness and False Detection Rates 


Before presenting an analysis of our detection statistics, we discuss the relationship between 
frequently used terms in the literature. After ADI processing of adaptive optics images, the residual 
intensity statistics is nearly Gaussian (Mawet et al. 2014, see also for detection statistics at 1-3 


X/D from the primary). A signal S embedded in Gaussian noise with standard deviation, cr =1, 
when detected using a signal threshold T, will be recovered at a completeness (or true positive rate) 
of C — G(S — T) = J^°_ s e~ x2 / 2 dx (with inverse function, G'(C) — S — T; Equation 1). This 

is because the real signal is redistributed due to the noise in a Gaussian manner (see Figure [4]). 
In other words, the completeness is the fraction of real (or injected) signal above a level, S', that 
is detected by a method. The FDR, which we measure as the number of bogus detections or the 
number of detections in the absence of signal, is theoretically FPR/(FPR+C), where the FPR is 
G(-T). 


Since there is a one-to-one correspondence between SNR improvement and contrast improve¬ 
ment, (Acontrast=2.51og 10 mag), we only plot contrast improvements. According to equa¬ 

tion 1, an improvement in SNR from 5 to 10 corresponds to improvement in completeness from 
50% to nearly 100% (to within 10 -6 ), for T— 5. Thus, completeness improvements quickly become 
insensitive to SNR improvements. Nevertheless, since it is interesting to show the improvement 
in the fraction of sources recovered, we compare the completeness for PC A and basic ADI at the 
threshold where MLOCI yields 50%. 

In Figures |6}fT0l we compare average detection statistics from our 10 NICI datasets for basic 
ADI, PCA and MLOCI for different reduction parameters. In Figure [6j we see that the perfor¬ 
mance of MLOCI is not very sentitive to P, the relative intensity of the negative sources, generally 
yielding contrast improvements of 0.3-0.4 mag. The plots are made for <a=2.15 and optimum pa¬ 
rameter settings W —16 and N— 2. However, for very large P (e.g. 10-100) we have found that the 
performance of MLOCI does deteriorate. Moreover, the counterpart PSFs begin to look similar to 
each other, instead of the science images that they are being matched to. This is because LOCI 
strongly biases the linear combination to match the negative sources in the counterpart images. 
The completeness improvement for MLOCI at the chosen threshold is usually 50%. The ADI and 
PCA performances are nearly equal to each other. 


From our reduction experiments, we can estimate the accuracy of the photometric measure¬ 
ments allowed by MLOCI. The systematic uncertainty of the photometry can be estimated by the 
change in the contrast before and after recovery, between the companion in the test sector to the 
two point sources in the reference sector. Since the reduction is done 10 times with the reference 
companions in different positions, we can also estimate the random uncertainty by calculating the 
standard deviation of the contrast after recovery between the two sets of companions. When an 
actual companion is detected in a real dataset, we can estimate photometric uncertainties and sys¬ 
tematic errors by injecting companions at other PAs at the same separation and performing MLOCI 
reductions to recover them. In the bottom panel of Figure [6j we see that while the photometric 
uncertainty is insentitive to P (variation of ^13% in these cases), the systematic error increases 
with increasing P. 

In Figures |6]-[8j we see that the contrast (and completeness) improvements are not very sensitive 
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to the width of the reference sector or the number of companions used as reference, or the separation 
of the recovered companions from the primary. However, Figure [T0| shows that for very small sky 
rotation (a =0.35) the contrast gain due to MLOCI is much larger (^1.4 mag). Here, PCA also 
shows a 0.5 mag gain over ADI. The variation of photometric uncertainties over these experiments 
is small. 

In Figure [TlJ we show attempts at measuring the source properties at the test locations after 
800 MLOCI reductions of the UY Pic data set when no actual signal has been injected. According 
to the previously defined detection criteria in §4, none of the measurements satisfied the detection 
crtieria and thus no false positives were detected. Moreover, the standard deviation of the SNR 
distribution was 1.0, showing that the false positive rate is not elevated above that expected. A 
higher false positive rate would indicate that while the reduction method may produce detections 
with seemingly high SNR, the noise is in fact not being calculated correctly. 

We also performed non-blind reduction experiments where the test companion and the reference 
companions are treated the same by MLOCI, using negative sources in all locations when creating 
the counterpart PSFs and/or positive sources when combining all the images. In these cases, 
the apparent contrast improvement was on average 0.75 mag. However, using the best algorithm 
parameters we were only able to reduce the false detection rate to 1%. Poor parameters could yield 
false detection rates > 20%. Thus, MLOCI must be kept blind about any suspected companion, 
except when informed by another blind process, which in our case is a basic ADI reduction. By 
the same logic, other versions of LOCI and PCA can also accidentally create false detections, 
since no restrictions on the allowed linear combinations are generally imposed when constructing 
reference images. This is an especially unpleasant possibility when trying to re-detect a companion 
marginally detected in an earlier epoch. We therefore strongly recommend that FDRs be provided 
for any reduction method. 


4.3. GPI Early Science Data 


Here, we describe the results of MLOCI and ADI processing of the GPI data on HD 135344 B. 
In Figure [3^, we show a de-rotated median combination of all the channels without PSF subtraction, 
but scaled in intensity by p 2 , and then bandpass filtered to isolate spatial features between 2-8 
pixels. We easily see the spiral arms, discovered in Muto et al. (2012), extending out to 0.9", and 
the cavity inside 0.2". However, we do not see their inner ring where the spiral arms begin. We do 
see a bright rim at 0.15" separation around the coronagraphic mask. Since, there is presumably 
no self-subtraction of the disk (no ADI) in this image, we use this image to identify reliable signal 
and background regions for MLOCI (Figure |3]l). 

In the regular ADI reduction (Figure^) we find some features at 110-170 mas separations 
with SNR=3-5 (AJ=8-9.5 mag), particularly two radially extended spots to the NE and SW. As 
discussed in Menard et al. (2015, in preparation), while these could prove to be accretion streams 
from the disk onto the star, this is currently uncertain because of the small amount of sky motion 
(1.7-2 xFWHM) at the relevant separation and similar PSF features in the unreduced images near 
mask (see Figure^). Moreover, GPI images of HR 4796 A (Perrin et al. 2014) also show similar 
residuals near the mask. 
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We find that particular caution needs to be exercised with MLOCI when there is signal in the 
FoV which is extended azimuthally, as in the case of HD 135344 B. When we try to target each 
spot separately with non-blind MLOCI, we obtain isolated point sources, admittedly with large 
FWHM, while the rest of the disk disappears (see Figures [3^). This demonstrates the ability of 
non-blind algorithms which allow free linear combinations of science images to produce misleading 
features, especially in the case where disk emission may exist. In the absence of disk emission, 
we have found that the most conservative non-blind methods create less false positives but still at 
elevated levels compared to ADI. 

However, the MLOCI algorithm for extended-sources does have the special ability to improve 
the SNR of any feature suspected in a dataset, be it mutiple companions or an entire disk. To 
demonstrate this potential of the algorithm, we use MLOCI (P= 3; extended-source reduction as 
described in §2 & §3) to optimize the SNR according to regions selected in Figure [3]l. Indeed, the 
spiral and the suspected streamers are much enhanced in the MLOCI reduction (see Figure^). 
When we flip the matching regions vertically so that they do not match with the real disk emission, 
MLOCI yields no disk (Figure [3]f). Thus MLOCI is not likely to yield false disk detections. 


5. CONCLUSIONS 


We have presented a version of the LOCI algorithm that significantly improves the SNR of 
both point and extended sources. This algorithm, called Matched LOCI or MLOCI, can be applied 
to both ADI and IFS data. The only difference being that, in the PSF construction step of image 
processing, the images have to be pupil aligned in the former and both pupil and spectrally aligned 
in the latter. The contrast gain over both regular ADI and PC A in NICI data was 0.3-1.4 mag. 
This is especially interesting, since it was shown that LOCI itself provides no gain over regular 
ADI processing for NICI data (Wahhaj et al. 2013a). While we found that LOCI could reduce the 
speckle noise in the final reduced image, it also reduced the signal of test companions. 


We have shown how MLOCI can be used in the detection of disks using GPI data, recovering 
new streamer-like features which need new observations before they can be confirmed. Thus MLOCI 
can be a powerful tool for the detection of companions and disks in the current direct imaging 
campaigns with GPI and SPHERE. 


While, for blind MLOCI the false detection rate is low, we found that for non-blind MLOCI 
for both NICI and GPI data, the false detection rate (FDR) can in some cases be much higher. 
The FDR is even more for when there is azimuthally distributed diffuse light in the FoV. Thus, in 
non-blind variants of LOCI and PCA algorithms high FDRs may also arise. Therefore, we strongly 
recommend that FDRs be measured whenever these methods are used. 


Our research has employed the 2MASS data products; NASA’s Astrophysical Data System; 
the SIMBAD database operated at CDS, Strasbourg, France. L.A.C. was supported by ALMA- 
CONICYT grant number 31120009 and CONICYT-FONDECYT grant number 1140109. L.A.C., 
M.R.S., H.C., and S.C. acknowledge support from the Millennium Science Initiative (Chilean Min¬ 
istry of Economy), through grant Nucleus P10-022-F. 
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Fig. 1.— The central steps of the MLOCI algortihm (see §2). Step 4, top row: For each science 
image a counterpart PSF is created. The PSF is such that its subtraction from the science preserves 
signal at 2 locations. Step 5, bottom row: Each science image and its counterpart PSF are de- 
rotated to align North Up. 
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Fig. 2.— The final steps of the MLOCI algorithm (see §2). Step 7, top row: Instead of the usual 
difference image combination in ADI, all science and PSF images are combined to recover injected 
sources in reference sector. Step 9, bottom row: Several such reductions created through step 7, are 
combined to create the final reduction. The final reduction preserves point source signal, generally. 
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Fig. 3. — All images are normalized by the RMS over the full images. They are oriented are North 
up and East to the left, (a) Simple stack of the HD 135344 B GPI J-band dataset scaled by p 2 
and spatially filtered to isolate features between 2-8 pixels in size, (b) Regular ADI reduction 
(c) Channel 35 (A=1.337/im) reduced with MLOCI X, enhancing only the SW spot, (d) MLOCI 
matching regions: Signal is preserved in the lightest regions. Noise is suppressed in medium blue 
regions. Darkest regions are not matched, (e) MLOCI reduction matched to regions in d. Region 
within 240 mas of center divided by 5, to improve stretch, (f) MLOCI reduction matched to map 
in d but flipped vertically. This demonstrates that it is difficult to create false disks by MLOCI. 
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Fig. 4.— Detection rates for the example case of 3cr signals and a detection threshold set at 2a. 
Here, we illustrate the relationship between False Positive Fraction (FPF), False Detection Rate 
(FDR), True Detection Rate (TDR) and True Positive Rate (TPR or Completeness) and how they 
depend on True and False Positives (TP shown in green and FP showin in red) and True and False 
Negatives (TN shown in grey and FN shown in yellow). The FDR, which determines the telescope 
time spent following-up bogus detections, should be minimized. This can be quite different in 
magnitude from the FPR. 
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Fig. 5.— Comparison of MLOCI reductions vs PCA and basic ADI reductions of NICI UY Pic 
data injected with simulated companions. The left panels in green show reductions for sky rotation 
equivalent to 0.35 xPSF FWHM (a =0.35) of motion at 0."5 separation. The blue panels to the 
right show the ones for 1.25xPSF FWHM of motion. The recovered companions can be seen at 
separations 0A5, OZ'6, OZ'7 and 0."8 with PA=360, 270, 180 and 90° respectively. The images have 
been normalized to the peak flux of the recovered companion at the top position (color scale is 
linear). The SNR improvement for MLOCI over PCA were 1.3 to 2.3 for a =0.25 and 1.3 to 1.9 
for a =1.25. The SNRs from PCA and basic ADI reductions were within 10% of each other. 
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Fig. 6.— Comparison of MLOCI to PCA and regular ADI for different pull factors (see § 2). We 
attempted to recover simulated companions with contrasts of ~9-ll mag at 0A5 separation. We 
compare the 95% completeness contrast and completeness improvements for regular ADI, PCA, 
and MLOCI reductions. The subtraction annulus width W is set to 16, while the artificial sky 
motion a was set to 2.15xFWHM. 
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Fig. 7.— Similar to Figure [6j a comparison of MLOCI to PCA and regular ADI for different 
subtraction region widths and <a=2.15 and P= 3. 
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Fig. 8.— Similar to Figure [6j a comparison of MLOCI to PCA and regular ADI for different 
separations and <a=2.15, W= 16 and P—3. 
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Fig. 9.— Similar to Figure[6j a comparison of MLOCI to PCA and regular ADI for different number 
of simulated companions and <a=2.15, W —16 and P— 3. 
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Fig. 10.— Similar to Figure [6j a comparison of MLOCI to PCA and regular ADI for different 
amounts of sky rotation and P =3 and W— 16. 
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Fig. 11.— The black dots show the measured signal to noise and PSF shape (same as fractional 
reduction in §4.1) of point sources injected at tests spots after 800 MLOCI reductions. The green 
dashed lines show the detection criteria defined in §4. The red dots show measurements for another 
800 reductions where no signal was injected. No false positives were found; the red dots lie outside 
the green box and reasonably far from the detection limits. Moreover, the standard deviations of 
the black and red SNR distributions are 1.2 and 1 respectively. 











